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We present an extensive study of angle-dependent transverse magnetoresistance in bismuth, with 
a magnetic field perpendicular to the applied electric current and rotating in three distinct crystallo¬ 
graphic planes. The observed angular oscillations are confronted with the expectations of semi-classic 
transport theory for a multi-valley system with anisotropic mobility and the agreement allows us to 
quantify the components of the mobility tensor for both electrons and holes. A quadratic temper¬ 
ature dependence is resolved. As Hartman argued long ago, this indicates that inelastic resistivity 
in bismuth is dominated by carrier-carrier scattering. At low temperature and high magnetic field, 
the threefold symmetry of the lattice is suddenly lost. Specifically, a 2n/3 rotation of magnetic field 
around the trigonal axis modifies the amplitude of the magneto-resistance below a field-dependent 
temperature. By following the evolution of this anomaly as a function of temperature and magnetic 
field, we mapped the boundary in the (field, temperature) plane separating two electronic states. 

In the less-symmetric state, confined to low temperature and high magnetic field, the three Dirac 
valleys cease to be rotationally invariant. We discuss the possible origins of this spontaneous valley 
polarization, including a valley-nematic scenario. 

I. INTRODUCTION 

Electric conduction in solids is affected by the appli¬ 
cation of magnetic field in a variety of ways. The most 
prominent is orbital magnetoresistance, which is the en¬ 
hancement of resistivity due to the Lorentz force suffered 
by charged carriers in presence of magnetic field. As early 
as 1928, Kapitza discovered that the electric resistivity of 
bismuth increases by many orders of magnitude upon the 
application of a large magnetic field [I]. The large orbital 
magnetoresistance is one manifestation of the extreme 
mobility of carriers in bismuth, itself a consequence of 
the lightness of electrons, the ultimate reason behind the 
singular role played by this elemental semi-metal in the 
history of scientific exploration of electrons in metals [S], 

During the last few years, bismuth has attracted new 
attention (A recent review can be found in ref.[5])- 
A number of intriguing observations on bismuth crys¬ 
tals exposed to strong magnetic fields have have been 
reported|3]-|S]. The angle-resolved Landau spectrum has 
been found to become exceptionally complex in high 
magnetic field. In spite of this complexity (and in con¬ 
trast to what was initially thought [U [5]), the spectrum 
resolved by experiment [5l [5 HTUI [HI is in agreement with 
theoretical expectations OH2] based on the band struc¬ 
ture of the system and its fine details. 

Open questions remain however. The three small 
pockets of Fermi surface residing at the L-point of the 
Brillouin zone host Dirac fermions with an extremely 
anisotropic mass becoming as small as one-thousandth 
of the free electron mass along the bisectrix axis. These 
three Dirac valleys are interchangeable upon a 27r/3 ro¬ 
tation around the trigonal axis. But, according to two 
sets of experimental studies [13J [14], the three Dirac val¬ 
leys become inequivalent at low temperature and high 
magnetic field. The origin of this spontaneous loss of 


threefold symmetry is yet to be understood. 

Metals hosting a small concentration of high-mobility 
carriers and displaying a large magnetoresistance have 
attracted much recent attention. In dilute metals such 
as WTe 2 [U! or Cd 3 As 2 [IS], resistivity enhances by many 
orders of magnitude upon the application of a magnetic 
field of 10 T. This is also the case of bismuth [B, 7} ITT] 
and graphite [T71 [18j, two well-known semi-metals. Both 
the amplitude of magnetoresistance and its held depen¬ 
dence have been put under scrutiny and are explored and 
discussed by experimentalists and theorists. 

In this paper, we present an extensive study of trans¬ 
verse angle-dependent magnetoresistance in bismuth and 
establish a detailed map of transverse magnetoresistance 
for all possible orientations of magnetic held from room 
temperature down to 2K and up to a magnetic held of 
14 T. The study aims to address two distinct issues. The 
hrst concerns the amplitude of magnetoresistance in a 
compensated semi-metal such as bismuth. Our results 
show that in any real material, the knowledge of all com¬ 
ponents of the mobility tensor is required to compute 
the magnitude of magnetoresistance. The emergence of 
a valley-polarized state is the second issue addressed in 
this study. 

Numerous studies of angular oscillations of magne¬ 
toresistance in strongly-correlated electron systems have 
been reported in the past. For three case studies, see 
Ref. [19l - l2T| . These investigations used angular magne¬ 
toresistance as a probe, knowing that it shows an ex¬ 
tremum when the peculiar topology of the Fermi sur¬ 
face modihes the dimensionality of the cyclotron orbit 
at a ’magic’ angle (dubbed either Lebed or Yamaji, af¬ 
ter those who conceived these commensurability effects). 
Our experimental configuration is different. The mag¬ 
netic field is rotated in a way to keep the magnetic field, 
B , and charge current, j, perpendicular to each other and 
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the transverse magnetoresistance is measured. In this 
configuration, the macroscopic Lorentz force between B 
and j, is constant. If the mobility were a scalar, no angu¬ 
lar variation would arise. Large angular oscillations are 
visible even at room temperature and in fields as small 
as 0.7T in bismuth|1 ,'1|. because carriers have drastically 
anisotropic mobilities. 

Our results provide an opportunity to test the semi¬ 
classic transport theory in a particularly constraining for¬ 
mat. Since the structure of the three electron ellipsoids 
and the single hole ellipsoid in bismuth are well-known, 
one can attempt to fit the experimental data by assum¬ 
ing reasonable values for the components of the mobility 
tensors of electrons and holes and establish their temper¬ 
ature and field-dependence. We find that, save for a num¬ 
ber of important details, our experimental results are in 
reasonable agreement with the expectations of the semi¬ 
classic theory. This paves the way to a quantitative un¬ 
derstanding of transverse magnetoresistance in bismuth 
with obvious implications for other semi-metals and di¬ 
lute metals. 

We also present a detailed study of the configuration 
in which the current is applied along the trigonal axis 
and the magnetic field is rotating in the (binary, bisec¬ 
trix) plane. We confirm the loss of threefold symmetry at 
low temperature and high magnetic field previously re¬ 
ported by transport [13] and thermodynamic studies |14j 
and find a boundary in the (field, temperature) plane sep¬ 
arating two electronic states. The underlying threefold 
symmetry of the zero-field crystal lattice is lost in the 
low-temperature-high-field state, but is kept in the high- 
temperature-low-field state. A phase transition between 
these two states is clearly detectable at low magnetic 
field. With increasing magnetic field, the transition shifts 
to higher temperature and becomes broader. We discuss 
possible origins of this phase transition and consider the 
available theoretical scenarios invoking Coulomb interac¬ 
tion among electrons [22] or lattice distortion induced by 
magnetic field reminiscent of Jahn-Teller effect [231. None 
of the currently available pictures provide an adequate 
description of the whole range of experimental facts. 


II. EXPERIMENTAL 

Most measurements were performed using a Quantum 
Design PPMS apparatus. These studies were comple¬ 
mented with measurements in Seoul using a two-axis 
home-made set-up in particular to check for any arti¬ 
fact resulting from misalignment. Three bismuth single 
crystals were cut in a cuboid shape with typical dimen¬ 
sions of (3 x 4 x 5 mm 3 ) and were used for mapping the 
magnetoresistance in three perpendicular planes. The 
typical residual resistivity of these crystals was 1 pfl cm. 
Taking a carrier concentration of n = 3 x 10 17 cm -3 for 
both hole-like and electron-like carriers and using the 
simple expression a 0 = ne{< p e > + < /q, >), the 
zero-field conductivity would imply that the sum of the 


average mobility of electrons and holes is as large as 
< p e > + < Hh >= 2 x 10 7 cm 2 V _1 s _1 . As we will 
see below, the order of magnitude is confirmed by our 
magnetoresistance data, but the mobility is very differ¬ 
ent for electrons and holes and along different orienta¬ 
tions. Samples of various cross sections (circular, square, 
triangular) were studied in order to check the effect of 
sample geometry on the loss of threefold symmetry. No 
significant difference between crystals of different origin 
(commercial vs. home-grown) was observed. 


III. ANGULAR OSCILLATIONS OF 
TRANSVERSE MAGNETORESISTANCE IN 
THREE PERPENDICULAR PLANES 

The Fermi surface of bismuth consists of a hole ellip¬ 
soid and three electron ellipsoids [24] . The hole ellipsoid 
has a longer axis, which is three times longer than the two 
other shorter axes and lies along the trigonal axis. The 
longer axis of each electron ellipsoid is about 14 times 
longer than the two shorter axes and lies in a (bisectrix, 
trigonal) plane, slightly tilted off the bisectrix axis. In 
the absence of magnetic field, charge conductivity in bis¬ 
muth is almost isotropic. But, behind this quasi-isotropy 
hides an intricate structure of opposite and compensat¬ 
ing anisotropies, which reveals itself by the application 
of a rotating magnetic field. 

Fig. 1 presents polar plots of transverse magnetoresis¬ 
tance as the electric current is applied along one crystal 
axis and the magnetic field is rotated in the plane per¬ 
pendicular to the applied current. The figure compares 
the data obtained for three perpendicular planes at a 
temperature of 30 K and a magnetic field of 0.5 T. The 
three polar plots illustrate how the structure of the Fermi 
surface leads to very different patterns in each case. 

When the current is applied along the trigonal axis 
(Fig. la) , magnetoresistance shows six fold angular os¬ 
cillations. The system has a C3 symmetry and remains 
invariant when the magnetic field rotates by 27r/3. In 
this configuration, the hole conductivity does not de¬ 
pend on the orientation of magnetic field. On the other 
hand, electrons show a large variation in their magne¬ 
toresistance. It becomes largest when the field is aligned 
along a bisectrix orientation, since for this orientation of 
magnetic field, Fermi velocity of electrons and the micro¬ 
scopic Lorentz force they feel are maximum. When the 
current is applied along the bisectrix axis and the field 
rotates in the perpendicular plane (Fig. lb), magnetore¬ 
sistance respects two angular symmetries. It is mirror 
symmetric (i.e. p(9 ) = p(—0) ) and it keeps the inversion 
symmetry(i.e. p(6) = p (tt + 0) ). When the current is 
applied along the binary axis and the field rotates in the 
third plane(Fig. lc), only the last [inversion] symmetry 
is kept. For the last two configuration, one expects a 
maximal magnetoresistance for holes when the magnetic 
field is along the trigonal axis. Added to the sum of the 
contributions of the three electron pockets, this leads to 
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T =30K B= 0.5T 


FIG. 1: Polar plots of transverse magnetoresistance as a function of the orientation of the magnetic field at T=30K and B=0.5 
K. The three plots show the data for three perpendicular planes, Left: held rotates in the (binary, bisectrix) plane; Middle: 
held rotates in the (binary, trigonal) plane and Right: held rotates in the (biscetrix, trigonal) plane. The electric current is 
always applied along the crystal axis perpendicular to the rotating plane. Upper panels show the projection of the Brillouin 
zone, as well as the hole and electron pockets of the Fermi surface in each rotating plane. 


complex patterns. As seen in Fig. lb and lc, the total 
magnetoresistance peaks at intermediate angles off the 
high-symmetry axes. 

These polar plots resemble those reported by Mase, 
von Molnar and Lawson half a century ago|25j. These 
authors presented a study of the galvanometric tensor of 
bismuth for an arbitrary orientation of the magnetic field 
for a single temperature (20.4 K) and a single magnetic 
field (0.576 T). The angular variation of our data is in 
good agreement with their data. However, the absolute 
amplitude of the magnetoresistance in our data at20 K 
and 0.5 T, is roughly two times lower. 

The extreme sensitivity of magnetoresistance in bis¬ 
muth upon the orientation of magnetic field permits one 
to use a rotating magnetic held to tune the contribu¬ 
tion of each electron pocket (i.e. each of three valleys) 
to the total conductivity in the configuration presented 
in Fig. la[TH]. This provides an interesting opportu¬ 
nity for “valleytronics” mm, an emerging held of re¬ 
search focused on manipulation of the valley degree of 
freedom. Such angle-dependent magnetoresistance is ex¬ 
pected to be observable in any multi-valley system with 
anisotropic valleys, as recently demonstrated in the case 
of SrMnBi2[5S]. What distinguishes bismuth, however, 


is the visibility of such angular oscillations at room tem¬ 
perature and in magnetic helds lower than 1 TP3j. This 
is a consequence of the large mobility of electrons, which 
exceeds 10 4 cm 2 V _1 s _1 at room temperature, and their 
very anisotropic mass. No solid other than bismuth is 
currently known to present such properties. 

The evolution of the angle-dependent magnetoresis¬ 
tance in bismuth when the field rotates in the (binary, 
bisectrix) plane was previous reported and the data wa 
found to be described by an empirical formula for multi- 
valley conductivityp3]. In this paper, extending the mea¬ 
surements to the two other planes, our aim is to see 
if transverse magnetoresistance in bismuth can be ex¬ 
plained for the whole solid angle with a single theoretical 
model based on a Boltzmann semi-classical picture. We 
will show that this is indeed the case and the empirical 
formula used in ref. m is an approximation of a more 
general formula for multi-valley systems. 

We performed an extended study on several crystals in 
a wide range of temperature (2K< T <300 K) and mag¬ 
netic field(B < 12 T). A selection of data for B = 0.5 T 
and T < 75 K are presented in Fig. 2-4. Fig. 2 shows 
the thermal evolution of the angular magnetoresistance 
when the current is applied along the trigonal axis and 
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FIG. 2: Transverse angle-dependent magnetoresistance for a 
current applied along the trigonal axis and a field rotating in 
the (binary, bisectrix) plane for different temperatures. The 
two panels present the same data plotted in a linear (top) and 
logarithmic (bottom) scale. 



FIG. 3: Transverse angle-dependent magnetoresistance for a 
current applied along the bisectrix axis and a field rotating in 
the (binary, trigonal) plane for different temperatures. The 
two panels present the same data plotted in a linear (top) and 
logarithmic (bottom) scale. 


the magnetic field rotates in the (binary, bisectrix) plane. 
This configuration is identical to the one studied in ref. 
m and the results are quite similar. The evolution of 
the angle-dependent transverse magnetoresistance with 
decreasing temperature for the two other planes of rota¬ 
tion is shown in Fig. 3 and Fig. 4. In all three cases, the 
evolution is smooth with increasing magnetoresistance as 
temperature decreases. This is a consequence of the fact 
that all components of the mobility tensor enhance with 
decreasing temperature. At zero magnetic field, the en¬ 
hancement in mobility with decreasing temperature leads 
to an enhancement of conductivity. In presence of a field 
as small as 0.5 T, on the other hand, it leads to an in¬ 
crease in resistivity with decreasing temperature. This is 
because the amplitude of the orbital magnetoresistance 
is set by the mobility and it exceeds by far the zero-field 
resistivity of the system in this “strong-field” limit. 

The thermal evolution of the morphology of the three 
curves is instructive. For the first configuration (J || trig¬ 
onal), there is no visible change in the structure of angle- 
dependent magnetoresistance with cooling. In the case 
of the other configurations, a qualitative evolution is ob¬ 


servable and new extrema emerge in p(9 ) curves as the 
system is cooled down. Note also the gradual saturation 
in amplitude of magnetoresistance below 7 K, indicat¬ 
ing that mobility has attained its maximum amplitude 
in this temperature range. 

For a quantitative treatment of the data, we need to 
revise the semi-classic transport theory applied to the 
case of multi-valley system with anisotropic valleys. 

IV. THE SEMI-CLASSIC TRANSPORT 
THEORY 

Following an earlier work by Abeles and Meiboom[29|. 
Aubrey [30 wrote the following equation for the charge 
conduction by a Fermi pocket with a carrier concentra¬ 
tion n in presence of electric and magnetic fields: 

j = a.E = nepE + p j x B (1) 

Here, e is the electric charge and a and p are the con¬ 
ductivity and the mobility tensors. In the absence of 
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FIG. 4: Transverse angle-dependent magnetoresistance for a 
current applied along the binary axis and a field rotating in 
the (bisectrix, trigonal) plane for different temperatures. The 
two panels present the same data plotted in a linear (top) and 
logarithmic (bottom) scale. 


magnetic field, Eq. 1 becomes the familiar expression: 
a(B = 0) = nefi. 

Abeles and Meiboom|29] argued that one can derive 
Eq. 1 starting from the linearized Boltzmann equation, 
which can be stated in the following manner (31) : 

+Urf k x B)^g k = - 9 ± ( 2 ) 

ae k n r 

Here, f k is the Fermi-Dirac distribution in equilibrium, 
g k represents the deviation from this equilibrium distri¬ 
bution, e k is the energey and vt the velocity of an electron 
with a wave-vector k and r is the relaxation time. In this 
picture, the charge current is defined as: 



Eq. 1 condenses all material-dependent parameters in 
the mobility tensor defined as: 


A = rm 


Here ?fi _1 is the inverse of the mass tensor. In the 
simplest case derivation r is a scalar. However, this is 
not necessary and one can easily generalize to a case with 
different relaxation times along x , y, z. 

Aubrey’s important contribution was to find that by 
defining a tensor B whose components are projections of 
magnetic held along the three perpendicular orientations, 
one can write a general solution to Eq. 1 in the following 
manner EDI: 


a = ne (A 1 + B) 1 (5) 

The components of the matrix B are: 

( 0 -b 3 b 2 \ 

B = I B 3 0 —B\ 

\-B 2 B 1 0 ) 

Here, Hi, B 2 and B 3 are the projections of the mag¬ 
netic held along the three principal axes: 



It is worthy to underline that no assumption has been 
made on the magnitude of the magnetic held. This is 
to be contrasted with treatments based on the Jones- 
Zener expansion, which are only valid in the weak-helcl 
limit (/ iB < l)[5ij. The Aubrey approach[501, on the 
other hand, is most appropriate in the strong-held limit 
(gB > 1 ). 

Now, to see the physics behind this picture, consider 
a spherical Fermi surface with an isotropic (i.e. scalar) 
mobility of g. In this case, the orientation of the mag¬ 
netic held has no importance. Let us assume it oriented 
along the z-axis (that is, B 3 = B 2 = 0 and B 3 = B ). In 
this case, the solution implied by Eq. 5 becomes: 

/ cr_L a H 0\ 

o = -<J H <J± 0 

V 0 0 °\\J 

Here, the transverse (crj_), Hall (an) and longitudinal 
(fj||) components of the conductivity tensor take the fol¬ 
lowing familiar expressions: 

ney neg 

= = (6) 

In other words, magnetic held modihes only the trans¬ 
verse conductivity and leaves the longitudinal conduc¬ 
tivity unchanged. Moreover, the Hall component out¬ 
weighs the transverse magnetoresistance in the high-held 
(/jB > 1) limit. Now, consider an ellipsoidal Fermi sur¬ 
face, with an anisotropic mobility: 

(lM 0 0 \ 

A = o ^2 o 

\0 0 g 3 ) 


( 4 ) 
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Rotating the magnetic field in the (1, 2) plane and 
measuring electric conductivity along the third axis, ac¬ 
cording to Eq. 2, one would expect to find: 

a = _ ne/j.3 _ 

33 1 + B 2 p, 3 (p 2 cos 2 9 + Hi sin 2 9) 

Here 9 is the angle between the magnetic field and the 
binary (i.e. 2nd) axis. When pi p 2 , this equation 
becomes identical to the empirical formula used to fit 
the data obtained for the field rotating in the (binary, 
bisectrix) plane [15]. 

Let us now turn to a multi-valley system where the 
total conductivity is the sum of the contributions by each 
valley. In the specific case of bismuth, one can write: 


CTtot = ^2 +a h . ( 8 ) 

i=l—3 

Here a h is the conductivity tensor of the hole pocket 
and of is the conductivity tensor of the electron pocket 
indexed i. 

The structure of electron ellipsoids in bismuth is com¬ 
plex. They have no circular cross sections and do not 
lie along any symmetry axis. Their mobility tensor, like 
their effective mass tensor I33j . has four distinct and finite 
components. In the case of the electron ellipsoid which 
has its longer axis in the ( x, z ) plane, this tensor can be 
expressed as: 


fm o o\ 

A? = 0 M 2 /m 

\ 0 p± p 3 ) 

The three electron pockets are equivalent to each other 
through a 27 t/ 3 rotation. Therefore the mobility tensor 
for the two other electron ellipsoids would be: 

A 2 = ^2ir/3 ' Al • B- 2 k/3 i A 3 = ^An/3 ' Al • B'iir/3 

Here, Rg is the rotation matrix for a rotation angle of 
9 around the trigonal axis : 

( cos 9 — sin 9 0\ 
sin# cos 9 0 I 

0 0 l) 

The hole pocket is an ellipsoid with a circular cross 
section perpendicular to the trigonal axis. It has identical 
projections in the binary and bisectrix planes. Therefore, 
there are only two distinct and finite components and its 
mobility tensor v can be written as: 

(vi 0 0 \ 

v = I 0 v\ 0 I . 

\0 0 v A ) 


Using this formalism, one can compute the components 
of total conductivity tensor for a given set of seven pa¬ 
rameters (n i_ 4 , ^13 and n). In order to compare with 
the experiment, one needs to invert the calculated tensor 
and obtain the relevant component of the resistivity ten¬ 
sor ( p = < 7 -1 ). In our particular case, we have measured 
resistivity, pn, p 22 , P 33 , along three principal axes. The 
link between pa and the components of the conductivity 
tensor can be written as: 


Pii — 


1 

(7a T 5 (7 a 


(9) 


Where: 


5(7 n = 


0’l20’23C r 31 + 0’i3(J2l0'32 — <722013031 — <733<712<721 
<722 <733 — 023(732 

( 10 ) 


5(7 22 — 


0’120’23( T 31 + < J i30'2i0'32 — 0’ll0’23( J 32 — 0"33CTi2(721 
<711(733 — (713(731 

( 11 ) 


5(733 ~ 


0’120’23( T 31 + < J i30'2i0'32 — 0’ll0’23( J 32 ~ 0 _ 22 0 'l30'31 
CT11CT22 — (712(721 

( 12 ) 


One may wonder if there is any room for harmless ap¬ 
proximations here. These expressions contain numerous 
off-diagonal components. Onsager relations imply that 
a ij(B) = —aji(—B). Agkoz and Saunders have detailed 
the restrictions imposed by Onsager reciprocity relations 
for all 32 crystallographic point groups [34j. They demon¬ 
strated that in the case of A7 structure of bismuth and 
antimony, the off-diagonal components of the conductiv¬ 
ity tensor can have both even and odd parts. In par¬ 
ticular, when the field is along the binary axis, there 
is a so-called ’Umkehr’ effect, leading to the inequality 
023 (-Bi) 7 ^ —<732(Bi). This means that the first two terms 
of the nominators in equations 10 to 12 do not cancel out. 

It is true that in a compensated semi-metal with an 
equal concentration of electrons and holes, the opposite 
Hall response of electrons and holes can cancel each other. 
In case of negligible off-diagonal components, one would 
simply find pa ~ . When the magnetic field is aligned 

along each of the three high-symmetry axes, one can ex¬ 
perimentally verify that p z j <C pa and thus a^ <C <7a- 
However, for arbitrary orientations of the magnetic field, 
the magnitude of the off-diagonal components is large 
enough to generate non-negligible effects. We found that 
in order to find a satisfactory fit, the full expressions of 
equation 9-12 are to be used. In particular, in the case of 
the third configuration, i.e. the field rotating in the (bi¬ 
sectrix, trigonal) plan], without taking 5au into account, 
even a qualitative agreement with the experimental data 
is impossible to obtain. 

In 1974, Siimengen, Tiiretken and Saunders[55] em¬ 
ployed this formalism to quantify the components of the 
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FIG. 5: Comparison between the data and the theoretical fit based on the optimal parameters at B=0.5 T and three different 
temperatures for : a) Magnetic field rotating in the (binary, bisectrix) plane; b) magnetic field rotating in the (trigonal, 
bisectrix) plane; c) magnetic field rotating in the (trigonal, binary) plane. 


mobility tensor in bismuth using the experimental data 
published a decade before by Mase and co-workers [2S] for 
a single temperature (20.4 K) and a single magnetic field 
(0.576 T). They found that experimentally-observed an¬ 
gular variation can be theoretically reproduced and the 
extracted components of the mobility tensor match those 
found by measuring components of resistivity, ptj in the 
low-field (pB -C 1) limit [551 HIT] . In the next section, we 
are going to use this procedure for our extensive set of 
data. 


V. THE MOBILITY TENSOR 

According to the theoretical frame described in the 
previous section, the angle-dependent transverse magne¬ 
toresistance can be fit using equations 9-12 with seven 
adjustable parameters. These are the four components 
of the mobility tensor of electrons, P\~a, the two com¬ 
ponents of the mobility tensor of holes, zq^ and the 
carrier density of holes, nh, which is equal to the sum 
of the carrier densities of the three electron pockets 
(n h = n e i + n e2 + n e3 ). 

We tried such fits and found that good fits with plau¬ 
sible parameters can be achieved. Polar plots comparing 
the experimental data and the expectations of the theory 


are shown in Fig. 5. This agreement establishes that the 
large anisotropic magnetoresistance of bismuth can be, 
mostly if not entirely, explained by the semi-classic the¬ 
ory. At lower temperatures and higher magnetic fields, 
the fits become less satisfactory. 

There are several reasons for the gradual inadequacy 
of the model as the temperature lowers and the magnetic 
field increases. First of all, Landau quantization becomes 
sharp enough to introduce additional angular structure 
not expected by the model used here. The second reason 
is the loss of the threefold symmetry, which will be dis¬ 
cussed in detail in the following section. A third is the 
apparent effect of the applied magnetic field on the com¬ 
ponents of the mobility tensor. This last feature appears 
only after putting under scrutiny the details of the fitting 
procedure. 

We found that in some cases, two different sets of fit¬ 
ting parameters could both yield satisfactory fits of com¬ 
parable quality. In order to minimize this uncertainty, we 
attempted to reduce the number of adjustable parame¬ 
ters by excluding the carrier density and the tilt angle as 
variable parameters. 

It has been known that carrier density in bismuth re¬ 
mains constant up to 50 K and then begins to increase 
as thermally excited carriers are introduced across the 
small gap at the L-ppint pH)] . We took a fixed value of 















































n e = rih = 3 x 10 17 cm -3 (the magnitude obtained from 
de Hass-van Alphen measurements [39]) for T < 50 K and 
let it evolve as a free parameter for temperatures above. 

The tilt angle, 9, between the longer axis of the electron 
ellipsoid and the crystalline bisectrix axis, generates the 
fourth mobility component, /i 4 . If the relaxation time 
tensor happens to be aligned along the high-symmetry 
crystals axes, the mobility tensor and the mass tensor 
should have identical tilt angles. It was measured to 
be 6.4 degrees by de-Haas-van-Alphen effect studies[39] 
and its controversial sign was settled by refined com¬ 
parison between the data on structural and electronic 
properties^. The recent study of angle-resolved Lan¬ 
dau spectruml9j found a tilt angle(~ 6.2° ) close to the 
previous reports. 

In the case of the mobility tensor, the tilt angle can be 
expressed as 35 36]: 

0 = \ arctan(——) (13) 

2 /i2 ^ 1^3 

Hartman found that the tilt angle of the mobility ten¬ 
sor is very close to this value and concluded that the 
relaxation time tensor has little or no tilt [3B]. We found 
that the best fits to our data point to a value close to 
6 .8° and in order to reduce the fitting uncertainty im¬ 
posed this as a constraint to our fitting procedure. 

Figure 6 and Figure 7 show the the temperature- 
dependence of the mobility components of electrons and 
holes extracted from our fits to the data at B = 0.5 
T. Our data is compared with zero-field values extracted 
from the galvanometric coefficients of bismuth in the low- 
magnetic field limit. Several remarks are in order. 

First of all, the sheer magnitude of the mobility of 
Dirac electrons in bismuth is remarkably large. At 10 
K, fii in bismuth becomes as large 1000 T -1 or 10 7 
cm 2 V _1 s _1 . This is slightly lower than what Hartman 
found at the same temperature. According to his results, 
fj ,i becomes as large as 10 8 cm 2 V _1 s _1 . Therefore, and 
this may be worth being recalled, electrons in bismuth 
are by far more mobile than carriers seen in any other 
three-dimensional solid |T5] . 

Second, the mobility tensor is anisotropic. Qualita¬ 
tively, this reflects the anisotropy of the band mass. 
The anisotropy is larger for electron, which have a more 
anisotropic mass, than holes. The mobility of elec¬ 
trons is lowest along the bisectrix, which is the direc¬ 
tion for which the electrons are heaviest. As for holes, 
they are lighter and more mobile perpendicular to the 
trigonal axis than parallel to it. However, the corre¬ 
spondence between the two tensors remains qualitative. 
The mass anisotropy of electrons, extracted from studies 
of quantum oscillations [9j [39! is as large as 200. The 
anisotropy of mobility does not exceed 40. This im¬ 
plies an anisotropic relaxation time tempering the mass 
anisotropy. 

It is instructive to compare our mobilities at B = 0.5 
T for various orientations of magnetic field to those re¬ 
ported by previous authors in the zero-field limit [3B, ‘171. 



FIG. 6 : Temperature dependence of the components of the 
electron mobility tensor, /ii, obtained by fitting our angle- 
dependent magnetoresistance data at B = 0.5T for a current 
applied along the trigonal axis. Also shown are the results re¬ 
ported for zero magnetic field by Hartman below 15 K[36] and 
by Michenaud and Issi above 77 K [37]. Solid lines represent 
a T ~ 2 temperature dependence (No data on p ,4 was reported 
in the latter case). 



FIG. 7: Temperature dependence of the components of 
the hole mobility tensor, Vi, obtained by fitting our angle- 
dependent magnetoresistance data at B = 0.5T for current 
applied along the binary or bisectrix axis. Also shown are 
the results reported for zero magnetic field by Hartman by 
Hartman below 15 K[3^ and by Michenaud and Issi above 77 
K|37) (No v 3 was reported in the latter case). 

As seen in Fig. 6 in the case of electrons, a reasonable 
match is observed between the different sets of results. 
In the case of holes, on the other hand, as seen in Fig. 
7, a difference is visible at low temperature: the mobil¬ 
ities extracted from data saturate to values significantly 
lower than what was found in the zero-field limit. This 
may be due to a difference in the ultimate carrier nrean- 
free-path in different samples. The residual resistivity in 
our samples is an order of magnitude larger than those 
studied by Hartman. It may also partially arise from a 
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FIG. 8: Temperature dependence of the components of the 
electron mobility tensor, /q, obtained by fitting our angle- 
dependent magnetoresistance data at B = 0.5 T for the three 
planes of rotation. The discrepancy increases with decreasing 
temperature. 


field-induced decrease in the magnitude of zq and zq. 

The effect of magnetic field on the various components 
of mobility tensor shows itself in another fashion. Fig. 8 
compares the magnitude of the components of the mo¬ 
bility tensor of electrons for the three rotating planes. 
As seen in the figure, the discrepancy remains within 
the error margin at high temperatures, but exceeds it in 
the low-temperature regime. Moreover, the discrepancy 
has a clear pattern. All /q _4 components become lower 
in the case of the third configuration, with the current 
applied along the bisectrix axis and the magnetic field 
rotating in the (binary, trigonal) plane. This indicates 
that the components of the mobility tensor are affected 
by the magnitude and orientation of magnetic field in the 
low-temperature limit. 

Hartman [36] noticed that the components of mobil¬ 
ity tensor present a quadratic decrease with temperature 
over a wide temperature window. This T~ 2 tempera¬ 
ture dependence of mobility (and consequently zero-field 
conductivity) is a hallmark of electron-electron scatter¬ 
ing. Long ago, Baber [5T] argued that in a metal with 
multiple electron reservoirs coupled to the lattice ther¬ 
mal bath, electron-electron scattering should give rise to 
a T 2 inelastic resistivity. In the case of metals hosting 
correlated electrons, this T 2 resistivity has been exten¬ 
sively documented and is one of the two main ingredients 
of the much-discussed Kadowaki-Woods ratio [12]. 

One aspect of our data remains currently beyond our 
understanding. In the first configuration, in which the 
current is applied along the trigonal axis and the field is 
rotating in the (binary, bisectrix) plane, the contribution 
of the hole pocket is expected to be independent of the 
orientation of magnetic field. Indeed, for this configura¬ 
tion, the hole conductivity becomes simply: 
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FIG. 9: Evolution of the angle-dependent magneto¬ 

conductivity with decreasing temperature. Each panel shows 
a polar plot of the inverse of magnetoresistance normalized 
to its maximum value at a given magnetic field. At high- 
temperature, the threefold symmetry of the underlying lattice 
is preserved. As the temperature decreases, this symmetry is 
lost above a threshold magnetic field, which decreases in am¬ 
plitude with cooling. 


_ 3nezq 

CT33 “ lT 


(14) 


Therefore, the fit to the experiment at a fixed mag¬ 
netic field B yields only a single parameter proportional 
to the right hand side of this equation and not an inde¬ 
pendent estimation of zq and zq. The mystery is that 
the zq and the zqzq found from the data obtained in this 
configuration are five times larger than what was found 
for the two other configurations. This discrepancy is well 
beyond our uncertainty margin. It remains the main fail¬ 
ure of the semi-classical model employed here and points 
to an additional transport process not taken into account 
in the picture employed here. 


VI. EMERGENCE OF A VALLEY-POLARIZED 
STATE 


Let us now turn our attention to the loss of valley sym¬ 
metry, a feature reported in two previous experimental 
studies,12 Hl|. In order to explore this feature, we have 
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performed extensive measurements with numerous sam¬ 
ples performed in two different laboratories (Paris and 
Seoul) to document the emergence of spontaneous valley 
polarization. 

When an electric current is applied along the trigo¬ 
nal axis, the three electron valleys remain degenerate. 
Of course, this degeneracy is lifted by a magnetic field 
applied along an arbitrary orientation in the (binary, bi¬ 
sectrix) plane. However, the C3 symmetry of the sys¬ 
tem implies that a 27 t/ 3 rotation of the magnetic field 
around the trigonal axis should not affect the physical 
properties, since the degeneracy remains lifted in exactly 
the same way. Therefore, a magnetic field rotating in 
the (binary, bisectrix) plane, should generate an angle- 
dependent magnetoresistance keeping a six-fold symme¬ 
try in a polar plot, resulting from the combination of 
the C3 symmetry and the inversion symmetry. This is 
indeed the case of our data at high temperature and/or 
low magnetic field. At low temperature and high mag¬ 
netic field, on the other hand, this symmetry was found 
to be lost in all the samples studied. 

Fig. 9 presents a detailed study on a bismuth single 
crystal with a square cross-section. Each panel of the 
figure shows polar plots of p^ 1 at a given temperature 
for different magnetic fields. In the following, p~ l would 
be called “conductivity” keeping in mind that, because 
of the contribution of the off-diagonal elements a ^ p _1 . 
At each temperature, the data is normalized to the max¬ 
imum value of p _1 to allow an easy examination of the 
evolution as the temperature decreases and the magnetic 
field increases. One can see that the sixfold symmetry 
clearly present at 100 K for all magnetic fields is lost 
at 5K for all magnetic fields. In the intermediate tem¬ 
perature range, the sixfold symmetry is present at low 
magnetic field but is lost at higher fields. 

The loss of threefold symmetry with cooling, originally 
reported in ref. EE has been confirmed in more than ten 
different crystals of bismuth with different shapes and 
cut from three different mother single crystals and us¬ 
ing different experimental set-ups in Paris and in Seoul. 
The low-temperature asymmetries found in various sam¬ 
ples were different from one another. In particular, the 
departure from threefold symmetry is more striking in 
samples with a square cross-section than in samples with 
a cylindrical or a triangular cross-section. On the other 
hand, the temperature and magnetic field thresholds for 
the loss of symmetry were roughly similar. 

Various possible experimental artifacts as the source 
of the loss of threefold symmetry can be ruled out. Two- 
axis rotation experiments performed have found a simi¬ 
lar loss of threefold symmetry [13:. Therefore, it cannot 
be a consequence of uncontrolled misalignment between 
the magnetic field, the crystal axes and applied electric 
current. The fact that the frequency of quantum os¬ 
cillations remains identical for different orientations of 
magnetic field [2] (see also below) rules out any possi¬ 
ble role of internal uncontrolled strain, since experiments 
under strain have shown that presence of strain would 




FIG. 10: Top: The ratio of p h ^ x to f (p^, + p b > 2 + p bi * 3 ) 
as function of temperature in three different single crystals 
of bismuth at B = 0.1 T. In all three cases, a sudden jump 
occurs at T ~9 K. Bottom: polar plots of p _1 for sample PI. 


have caused a detectable variation among the three elec¬ 
tron pockets @3- Our bismuth crystals are twinned and 
therefore in addition to the dominant domain, there are 
three minority domains. The presence of these minor¬ 
ity domains generates a secondary set of Landau peaks 
visible by the Nernst measurements [TO]. But, we cannot 
think of anyway their presence can distort the threefold 
symmetry of the dominant domain, the only source of 
magnetoresistance in our study. 
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FIG. 11: Polar plots of resistivity for sample PI at different 
temperatures. Each panel shows a different magnetic field. 
An additional sub-structure arises by occupation and evacu¬ 
ation of Landau levels with the rotating magnetic field. The 
angular pattern becomes simple again at 8 T as no more Lan¬ 
dau level is occupied or evacuated during the rotation. 
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FIG. 12: The evolution of the jump observed in the ratio of 

PiTini to kiPbiil + PbLi + Pwns) ( see Fi S- 10 ) with increasing 
magnetic field in sample PI using the data plotted in Fig. 11 . 
As the magnetic field increases, the transition becomes wider 
and shifts to higher temperatures. 


By following the evolution of this loss of symmetry, we 
succeeded to identify a visible signature of a phase tran¬ 
sition associated with this loss of symmetry at low mag¬ 
netic fields. It is in the low-held window that an abrupt 
change becomes visible. As seen in Fig. 10, at B = 0.1 
T, the magnitude of magnetoresistance for magnetic held 
aligned along one bisectrix axis suddenly deviates from 
its average value at a temperature. Because of the resid¬ 
ual experimental misalignment, the amplitude of resistiv¬ 
ity for magnetic held oriented along each of the three bi¬ 
sectrix axes is slightly different. However, as seen in the 
hgure, this departure from unity is almost constant at 
higher temperature. Suddenly, below a held-dependent 
threshold temperature, an abrupt deviation on top of the 
experimental misalignment can be clearly detected. As 
seen in the hgure, in two other samples studied at B = 0.1 
T, a similar abrupt jump occurring at an almost identical 
temperature was found. This implies that the boundary 
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marking the loss of threefold symmetry is an intrinsic 
property of bismuth. 

In two samples PI, a cylinder with a circular cross 
section, and P2, a cuboid with a cleaved square cross 
section, we made extensive measurements up to 8 T in 
order to follow the evolution of this jump. The data for 
sample PI is shown in Fig. 11. As the magnetic field 
rotates, Landau levels are filled and evacuated and this 
adds additional substructure to the background. This 
substructure evolves as the magnetic field increases and 
finally vanishes at 8 T since no more filling or evacuation 
occurs. Nevertheless, the jump associated with the loss 
of threefold symmetry can be extracted by the same pro¬ 
cedure. The result is shown in Fig. 12. There are two 
distinct visible effects of the magnetic field on the phase 
transition. With increasing magnetic field, the transition 
shifts to higher temperatures and becomes wider. In a 
field as high as 8 T, the loss of symmetry occurs over a 
temperature range as wide as 20 K. On the other hand, 
the magnitude of the maximum deviation from unity does 
not change with magnetic field and remains a few percent 
in the whole range of investigation. 

Using this data, we can construct a phase diagram 
drawing the boundaries of the state in which the sym¬ 
metry of the three valleys is lost. This phase diagram 
is shown in Fig. 13. Each symbol represents the posi¬ 
tion of the intersection between the low-temperature and 
high-temperature behaviors detected in Fig. 12. The C3 
symmetry is preserved in the low-field-high-temperature 
region and lost in the high-field-low-temperature region. 
A very similar phase diagram is found for sample P2 
which is a cuboid with a square cross section. As seen in 
the figure, while the cross section does have a signature 
in the ultimate pattern of angular magnetoresistance in 
the low-temperature-high-field state, it does not modify 
the boundary in a significant way. This provides further 
evidence that this transition is intrinsic to bulk bismuth. 
We did not find any evidence for hysteresis in these mea¬ 
surements. 

We also studied the quantum oscillations of resistivity 
(the Shubnikov- de Haas effect) for different orientations 
of magnetic field along the three binary and the three 
bisectrix axes . We found that while the frequency of 
oscillations is identical for the three equivalent orienta¬ 
tions, there is a detectable difference in the amplitude of 
oscillations. This result is similar to what was reported 
by Kuchler et al. in the case of the quantum oscillations 
of magnetostriction ll4| and the quantum oscillations of 
the Nernst coefficient [8] . While there is no visible dif¬ 
ference between the angle-resolved Landau spectrum of 
the three electron valleys, the intensity of the signal gen¬ 
erated by emptying a Landau level is different among 
valleys. As seen in Fig. 14, the amplitude of oscilla¬ 
tions obtained after subtracting a smooth background 
are different for three nominally equivalent orientations 
of magnetic field. However, the oscillations have identical 
frequencies. Moreover, there is no visible difference in the 
oscillations of the second-derivative of resistivity, which 




FIG. 13: Top: Comparison of the angular variation of normal¬ 
ized conductance in two samples with different cross sections. 
In the high-temperature-low-field state, they both display 
the threefold symmetry, but in the low-temperature-high-field 
state they loose it with different patterns emerging in the 
valley-polarized state. Bottom: Variation of the critical tem¬ 
perature for the loss of threefold symmetry with magnetic field 
in the two samples. Vertical (horizontal) error bars indicate 
the width of transition as the temperature(magnetic field) is 
swept at a fixed magnetic field (temperature). The boundary 
between the ordered (less symmetric) and disordered (more 
symmetric) states in the (field, temperature) plane is similar. 


wipes out any field-linear magnetoresistance. Therefore, 
it is tempting to conclude that the difference in amplitude 
of quantum oscillations for the three valleys is roughly 
linear in magnetic field, like the monotonous background. 

To sum up, this study confirms the emergence of valley 
polarization in low temperature and high-magnetic field 
reported previously pT51. . In addition and for the first 

time, an abrupt change reminiscent of a phase transition 
can be clearly seen in the data and the boundary of the 
valley-polarized state (for a magnetic field perpendicu¬ 
lar to the trigonal axis) can be established by following 
the evolution of this anomaly with field and tempera¬ 
ture. Thus, there is now a body of experimental results 
pointing to valley polarization in bismuth, a phenomenon 
not still understood. While several theoretical ideas have 
emerged, none of them offers a satisfactory explanation 
of the whole spectrum of the observations. 

One appealing line of thought is to follow the idea 
of valley nematicity proposed for quantum Hall sys¬ 
tems with multiple eccentric valleys [22]. Abanin and 
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co-workers argued that in quantum Hall systems, i.e. 
two-dimensional systems subject to quantizing magnetic 
fields, if the valleys present an anisotropic effective mass, 
unequal occupation of valleys would save exchange en¬ 
ergy among electrons. As a consequence, a ground state 
in which one of the valleys is preferably occupied by elec¬ 
trons, dubbed a valley-nematic state, can be induced by 
magnetic field. In this scenario, the Coulomb interac¬ 
tion and the mass anisotropy in each valley play an im¬ 
portant role. The scenario is conceived mainly with a 
two-dimensional system with C4 symmetry such as AlAs 
heterostructures [33]. Bulk bismuth has a C3 symmetry 
and is three-dimensional. On the other hand, it has a 
mass anisotropy forty times larger than AlAs, making it 
a potential candidate for this scenario. 

Experiment, however, indicates that the three 
electron-like pockets do not differ in their size. The iden¬ 
tical frequency of quantum oscillations implies identical 
valley occupation. This is not conform to the expecta¬ 
tions along this line of thought. The three valleys may 
differ in the density of states near the chemical potential. 
This is what governs the magnitude of quantum oscilla¬ 
tions as well as the number of electrons participating in 
charge transport. But, how can the three valleys remain 
identical in size, but be different in their density-of-states 
near the chemical potential? There is no satisfactory an¬ 
swer to this question. However, Krichler et al. m recall 
the case of disordered semiconductors, in which the com¬ 
bination of Coulomb interaction and disorder open a gap 
in the immediate vicinity of the chemical potentialj45j. 
It is still unclear how one can use this idea in the case of 
a multi-valley metal such as bismuth. 

Another possibility is a field-induced lattice distortion 
driven by electron-phonon coupling. This idea, proposed 
by Mikitik and Sharlai[23], is reminiscent of Jahn-Teller 
distortion in insulators [IQ. According to the Jahn-Teller 
theorem, a number of exceptions notwithstanding, a geo¬ 
metrical configuration of atoms becomes unstable in pres¬ 
ence of electron degeneracy[47], Lowering symmetry al¬ 
lows the molecule to get rid of the costly electronic de¬ 
generacy. In the case of bismuth, the degeneracy of the 
three valleys has an energy cost, which may be reduced 
thanks to lattice distortion. There is currently no evi¬ 
dence for such a field-induced lattice distortion. There 
has also been no attempt to pin it down if it exists. 

Can a phase transition with a critical temperature of 
9 K be generated by a magnetic field as small as 0.1 T? 
To address this question, one should recall that the cy¬ 
clotron energy is inversely proportional to the mass of 
electrons. While, for free electrons, the cyclotron en¬ 
ergy at 0.1 T is only 0.13 K, in bismuth, where electrons 
are one thousand times lighter along the bisectrix, the 
cyclotron energy becomes orders of magnitude larger at 
the same magnetic field and a small magnetic field can 
cause a relatively robust instability. 
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FIG. 14: Top: Comparison of quantum oscillations of magne¬ 
toresistance (the Shubnikov-de Haas effect) for magnetic field 
along three equivalent bisectrix axes. The amplitude of the 
oscillating signal, obtained after subtracting a smooth back¬ 
ground for each angle is different along the three orientations. 
The difference becomes undetectable in the second derivative 
of magnetoresistance along the three different orientations. 
Bottom: The angle-resolved Landau spectrum revealed in a 
color plot of d 2 p/dB 2 keeps the threefold symmetry. 


VII. CONCLUDING REMARKS 

Decades ago, Wolff wrote: “The metal bismuth is 
one of the most peculiar and intriguing of solids” [33] . 
This peculiarity shows itself in numerous physical prop¬ 
erties ranging from its large diamagnetism to its remark¬ 
able thermoelectric figure of merit. In the case of our 
study, three basic facts are relevant: i) Carriers in bis¬ 
muth are extremely mobile; ii) This mobility is extremely 
anisotropic; iii) The anisotropy is inverted for electron¬ 
like and hole-like carriers. There are other solids display¬ 
ing these features, but the magnitude of mobility as well 
as the intricate complexity of its anisotropy puts bismuth 
in a league of its own. 

The first part of this study documents the evolution 
of angle-dependent magnetoresistance with temperature 
and magnetic field and finds that the complex angular 
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variation is within the reach of a semi-classical treatment 
if the mobility is taken as a tensor. This paves the way 
for a quantitative understanding of magnetoresistance in 
bismuth at arbitrarily large fields. Upon the application 
of a magnetic field as large as 65 T along the trigonal axis, 
bismuth remains metallic with an electric resistance in¬ 
creasing by six orders of magnitude^. This million-fold 
increase is even larger than what has been recently re¬ 
ported in \\ I'e/TT) in a comparable field range. How¬ 
ever, in contrast to the latter and what is expected for 
a simple compensated metalpI5]. the magnetoresistance 
does not present a quadratic field-dependence and cannot 
be described by a simple power law over an extended field 
window. One identified reason for this deviation is the 
field-induced change in carrier density at high magnetic 
fields. When the magnetic field exceeds 0.5 T, the de¬ 
generacy of Landau bands leads to a modification of the 
carrier number, while preserving the charge conservation. 
This feature, which has visible signatures in the high- 
field Landau spectrum [5] , should be taken into account in 
any realistic description of high-field magnetoresistance. 
Quantifying the components of mobility of electrons and 
holes in the low-field limit is a beginning. Understand¬ 
ing magnetoresistance at arbitrarily large magnetic field 
cannot be achieved by a toy model of a compensated 
semi-metal with a scalar mobility. Future studies will 
tell if one can describe the field dependence of magne¬ 
toresistance with realistic assumptions on the amplitude 
of the mobility components and their field dependence. 

The second part of this study documents the emer¬ 
gence of a valley-polarized state at low temperature and 
high magnetic field. Currently, this state, in which the 


electron fluid loses the threefold rotational symmetry of 
the underlying lattice, is poorly understood. Nematic- 
ity, a subject of intense attention in strongly-correlated 
electron systems, is associated with the loss of rotational 
symmetry by the electron fluid |50j. It is tempting, and 
maybe even legitimate, to qualify the valley-polarized 
state found in bismuth as valley-nematic. However, one 
shall not forget that the picture drawn by experimental 
results does not easily fit in the most straightforward de¬ 
scription of valley nematicity. In the ordered state, the 
three electron valleys remain identical in size but differ 
in their density of states at the Fermi level. There is no 
satisfactory understanding of how this can be achieved. 

To sum up, elemental bismuth hosts extremely mobile 
and very anisotropic carriers, which scatter off each other 
in a complex, yet increasingly documented, manner. This 
looks like an appealing platform to explore the possible 
existence of electronic equivalents of liquid crystals. 
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